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O ■ Abstract 

In the present paper we investigate the Hawking radiation of asymptotically non- 
flat dyonic black holes in 4D Einstein-Maxwell-dilaton gravity in semi-classical ap- 
proximation. We show that the problem allows an exact analytical treatment and 
CN ! we compute exactly the semi-classical radiation spectrum of both non-extremal and 

extremal black holes under consideration. In the high frequency regime we find that 
the Hawking temperature does not agree with the surface gravity when the magnetic 
charge is nonzero. Even more surprisingly the Hawking temperature is independent of 
the black hole intrinsic characteristics, as the mass and magnetic charge, and depends 



only on the linear dilaton background parameter. 
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1 Introduction 

The Hawking radiation of black holes is an emblematic effect in the quantum field theory in 
curved spacetime pQ, [2]. This effect, lying on the wedge of classical and quantum gravity, 
reveals the deep connection between the black hole physics and thermodynamics. In the 
semi-classical approximation, the spectrum of the Hawking radiation can be obtained by 
computing the Bogoliubov coefficients in two different vacua and matching them appropri- 
ately [2 J . Another procedure which, in asymptotically flat spacetimes, gives the same result 
as the Bogoliubov coefficients method, is to compute the absorption and the transmission 
(or reflection) coefficients of waves defined at the asymptotic regions [3]. 

Ever since its discovery the Hawking radiation continued to be a hot area of research. The 
importance of the Hawking effect for the fundamental physics stimulates its investigation 
for various black hole solutions with different structures and asymptotics. Unfortunately, 
the wave equations in black hole spacetimes can not be solved analytically and this makes 
the full study of the Hawking radiation hard. Only in special cases we are able to solve the 
wave equations exactly. The cases when the wave equation is exactly solvable are important 
because they enable us to study the Hawking radiation in detail and, for example, allow us 
to compute the radiation spectrum exactly. In the present paper we consider one such case 
when the exact analytical treatment of the Hawking radiation is possible. More precisely 
we study the Hawking radiation of a class of asymptotically non-flat dyonic black holes in 
the 4D Einstein-Maxwell-dilaton gravity. 

The Einstein-Maxwell-dilaton gravity in 4-dimensional spacetime is described by the 
following action 

S = j^ J Av^ [K - 2<rv MV >V„v> - e- 2 ^F, u F^} , (I) 

where TZ is the scalar curvature with respect to the spacetime metric g^ u , F^ u is the electro- 
magnetic field, and tp is the scalar dilaton field with a coupling constant a. In the present 
paper we are interested in the black hole solutions of ([I]). More precisely we consider the 
following dyonic black hole solution for a = 1 found in [I]: 

■>~ 2 ( '" ''- ){ ' ~ r+) dt 2 + - ^ -dr 2 + r r(d6 2 + sin 2 6d<J) 2 ) } (2) 

(r — r_)(r — r+) 

(3) 

(4) 

(5) 

where r_, r + and tq are constants. Here $ e and $ m are the electric and magnetic potentials 
and the Maxwell 2-form is given by 
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F = d$ e A dt + e 2v * (d® m A dt) (6) 

with * being the Hodge dual. 

The solution describes an asymptotically non-flat dyonic black hole with inner and outer 
horizons at r = r_ and r + , respectively. The electric and magnetic charges are given by 



In the limit P — > (r_ — > 0) we recover the pure electrical linear dilaton black hole 
solution [5] (see also [6]). It should be noted that the parameter r (or equivalently the 
electric charge Q) is associated not with a specific black hole, but rather with a given 
linear dilaton background. The linear dilaton background solution is obtained by setting 
r + = r_ = in our solution. 

In order to find the physical mass of the black hole solution under consideration we 
use the quasilocal formalism [TJ. Since our spacetime is not asymptotically flat, a suitable 
substraction procedure is needed to obtain a finite mass. In our case the linear dilaton 
background is the most natural and unique choice for the substraction background. The 
explicit calculations give the following result for the mass jl] 

M=^(r + + r_). (8) 

The surface gravity is given by 

k = r -±^. (9) 

In the case with r + = r_ we obtain an extremal black hole solution with zero surface 
gravity. 

Following j8] one can formally derive the first law for the black holes under considerations, 
namely 

dM = ^d (If) + $y P (10) 

where A^ = 47rror + is the horizon area and $^ = * /^jr- is the magnetic potential evaluated 

on the horizon r = r+. Here the parameter r$ (respectively Q) is kept fixed since it is 
associated with the background. General comments about the black hole thermodynamics 
in asymptotically non-flat spacetimes will be given in the last section. 

In the next section we study the Hawking radiation of the asymptotically non-flat black 
hole solutions presented above in semi-classical approximation in both the non-extremal 
and extremal case. The Hawking radiation of pure electrical linear dilaton black holes (i.e. 
corresponding to P = 0) was studied in semi-classical approximation in [9]. 



2 Hawking radiation in semi-classical approximation 

2.1 Hawking radiation of non-extremal black holes 

In order to study the Hawking radiation of our black holes we consider a test scalar field ip 
satisfying the wave equation 



□V> = o, 



(ii) 



where □ is the curved spacetime D'alambert operator. 

For the static spherically-symmetric metric 02]), the D'alambert operator takes the fol- 
lowing explicit form: 



□ 



d^yFgg^d, ) 



TqT 



1 «r^ - v- _. ^ i , j -. v 



rd tt + — d r [(r-rJ)(r-r + )d r } + — ^ (e>(p) 

(r— r-)(r — r + ) r^r r$r 



where A(o,tp) is the Laplace operator on the unit sphere S 2 =1 . So the scalar wave equation 
dip = multiplied by r r becomes: 



r 2 r 2 



-4V + d r {(r-rJ)(r-r + )d r \ ij) + A (e>(p) ip = . 



(12) 



(r— r_)(r — r + ) 

Consider now a harmonic eigenmode as a partial solution of (fl2l) in separate variables: 

il> ulm (t, r, 9, <p) := R ul (r)Y lm (9, <p)e~ luJt (13) 

where Yi m (8, cp) are the spherical harmonics. Then the Fourier coefficients R w i(r) satisfy the 
equation: 



d_ 
dr . 



(r-r-)(r-r + )Rui(r) 
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(r—r-)(r—r+) 



1(1+1) 



iO(r) = 



(14) 



where Go : = r u is a dimensionless frequency. 
After the substitutions: 



z := 



r + — r 



r + — r- 
the equation for Z u i(z) := R u i(r) becomes: 

z(l-z)Z ul (z) + (l-2z)Z wl (z) + 



z :-" 



r_i_ — r_ 



> 1 



Go 2 (z-z ) 2 
z(l-z) 



+ 1(1+1) 



Zui(z) 



0. 



(15) 



(16) 



An appropriate substitution like Z u i(z) := z p (z — l) q h UJ i(z) will help us to obtain a more 
familiar linear differential equation: 



z(l-z)h u i(z) + [l + 2p - 2(p+q+l)z} h w i(z) 
(p+q+l)(p+q) + Gj 2 + 1(1 + 1) - 



., uf m p 2 +z%u 2 q 2 +(z -l) 2 u 2 
z \ — z 



h u i(z) = 



where p and q are determined by the conditions which eliminate the simple rational 
fraction 1/z and 1/(1 — z) in the coefficient in front of h w i(z). There are four combinations 
for p and q, that can do it, but the most convenient one is to take p := izqCj and 
q := -i(z -l)Q\ 

z(l-z)Mz) + [l+2iz CJ - 2(l+iw)z] Mz) - [iw - 1(1+1)] M*0 = . (17) 

This is just the hypergeometric equation and after identifying ffTTj) with the canonical form: 

z{l-z)h u i{z) + [c-{l+a+b)z]K A {z) -abh u i(z) = 0, (18) 

one can easily obtain the canonical parameters: 

a= 1/2 + i(uj+\ ajl ) 

6=1/2 + i{CJ-\ wl ) , where A^ := ^u 2 - (Z + l/2) 2 . (19) 

c = l + 2iz u) 

In our dimensionless variable (fT5l) the singular point z = 1 corresponds to the inner 
horizon r = r_ while the singular point z = corresponds to the outer horizon r = r + . We 
are interested in the general solution around the outer horizon r = r + that can be continued 
to the spatial infinity z — > — oo (r — > +oo): 

Ki{z) =C X F (a, b, c; z) + C 2 z 1_c F (1+o-c, 1+6-c, 2-c; z) . (20) 

Consequently for 

Z ul (z) := z^iz-iy^-^h^z) , (21) 

we obtain the expression: 
Z ul (z) = (1 - z)-^"- 1 ^ [d^z^F (1/2 + 1(0;+^,), 1/2 + zCw-Ao,,), l+2iz u- z) + 

+ C 2 (-z)~^F (1/2 - i [(2z -l)u) -\J ,1/2-2 [(2z -l)cu +X U}1 ) , l-2*z u;; z)] . 

using (1T91 and (12TJJ). 

For the asymptotic of (T20]) near the outer horizon (i.e. z —¥ 0_) the zero order expansion 
of F (a, b, c; z) is enough (F(a, b, c, 0) = 1). Also lim z ^ _(l — z) -4 ^ ^ 1 ^ = 1. Then we have 

Z ujl (z^0 + ) ~ Ci(-zy zo * + C 2 (-z)- iZ0 * . (22) 

For physical interpretation it is appropriate to define a new real spatial variable x, by the 

relation: 

r ~ r+ x 

exp (23) 



r+ - r_ z r 

with r — > r + =>• x — >• —00 and r — >• +00 =>■ cc — >• +00. Now taking into account that x = x(r) 
is an uniformly growing function and replacing (T231 in (j2"2l . the asymptotic solution (|2"2"1) 
multiplied by e~ lw4 can be considered as a superposition of an out-going and an in-going 

wave 



Ru,i(r -¥ r , )e 



-iuit 



A out e l ^ x - t] + A in e~ iulix+t) , where 



A ou t '■— C\ 

Ai n := C*2. 



In the case I = 0, Y 00 (9,(p) = 1, and the upper expression is just the solution for the 
eigenmode 4>uioo(t, r , Q, f) for r = r + . 

At the spatial infinity (r —¥ +00 , 1/z —¥ 0), the asymptotic solution can be written by 
using the known relation between hypergeometric functions F(a, b, c; z) and F(al ', b', c'; 1/z). 
Taking only the leading order expansion of F(a', b', c'; 1/z) with respect to 1/z, we have the 
following expressions: 



F(b, c;z — >■ —00) 



r(c)r(6-q) 

T(b)T(c-a) 



x-q 1 r(c)r( a -6) b 

' + T(a)T(c-b) [ ] ' 



-zY' c F (l + a-c,l + b-c,2-c; z ^ -00) ~ 



(24) 



T(2-c)T(b-a) 



-z)~ a + 



r(2-c)r(o-6) 



T(l+b-c)T(l-a) y ~' ' r(l + a-c)r(l-6) v ' ' 
Applying the transformations (12"1"|) to (12TJJ) . replacing h w i{z) in (12T1) and finally taking 
into account that lim [— z/(l — z)\ = 1 , the general asymptotic solution becomes: 



1 



Z UI (Z -»- -OO) ~ -= [S^-*)**"' + SinC-^)-^'] 



— z 



(25) 



where 

r(c)r(q-6) r(2-c)r(a-6) 

-Dout : — ^1TT7 wT7 7T-r<-^23rT- tttt" — , £> 



) J -'in • 



r(o)r(c-6) T(l+a-c)T(l-b) 
Following ( 123|) and applying the substitution 

A^z 



c r(c)T(b-a) ]C ^ r(2-c)r(6-o) 



fe 



U)/ 



^O 



T(6)r(c-o) r(i+6-c)r(i-o) 



(26) 



when A^z G M, the radial function in x- variable can be considered again like a superposition 
of ID modes with wave vectors ±A; W ; ((— z) ul = e ±tkbjlX ): 



Ruji(r -¥ +00) 



r _ r _ 



[B out e^ x + B m e-^ x ] 



(27) 
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T(-2iX 



Ull, 





r(i + 2^ w) A out 
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_r(l/2 + i(d;H-A^))r(l/2 + z[(2^- 

T(l — 2izou) 


-l)w 
4- 


+■ A^]) 


1 r(i/2 


-i{u-\ wl ))T(l/2- 
T(l + 2£z w) A^ 


J [(20o 


-l)w- 
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_r(i/2 + z(w- 


-A^))r(l/2H-i[(2^ - 
T(l — 2iz u) 


-l)w- 

Ain 


-\,,]) 



(28) 



T(l/2 - z(w + A^)) r(l/2-i [(2z -l)cD + A^]) 
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Black hole radiation is considered as a specific boundary condition when only an out- 
going mode at the spatial infinity exists, B in = 0. This condition determines the ratio of 
the coefficients Ai n /A out or the reflection coefficient R: 



R 



\Ai. 



\A, -I 2 



'■out | 



r(l+2^ w)| 2 |r(l/2 - *(w+A^))| 2 |r(l/2 - i[(2z -l)u + A 



ul\ 



Bin—0 



|r(l-2i^)| 2 |r(l/2 + *(w-A^))| 2 |r(l/2 + i[(2z -l)u - A w/ ] 



Complex conjugation and the Euler's reflection formula for the Gamma function give us 
the final result for the reflection coefficient on the outer horizon: 



R 



In the special case r_ 



cosh(7r(o;-A^))cosh(7r((2^o-l)w-A a; i)) 
cosh(7r(d; + A^)) cosh(7r((2z -l)w+A^)) ' 

- =>- z = 1 we recover the result of [9] 
cosh 2 (7r(d)-A a; /)) 



R 



cosh 2 (7r(a; + A aj /)) 



(29) 



(30) 



For high frequencies Co ~^> 1 + 1/2 (=>- X^i ~ to) and u ^> -^r we obtain 



N : 



R 



cosh(27ru;) cosh(27T2;ow) 



1 - R cosh(2n(z — l)uj) 

i l + exp(-47r(z + l)w) 



cosh(27r(;a:o + l)w) 
cosh(27r(,2 — l)u) 



-i -l 



^Aituj 



1 + exp(— 4tt(zo — 1)oj) 



We identify the Hawking temperature from N ps e 



-AlTUJ 



-&1TLJ 



e t h which gives 



(31) 



T, 



H 



4nr 



(32) 



As one can see the black hole temperature derived in semiclassical approximation does 
not agree with the surface gravity, i.e. Th ^ ^ where the surface gravity k is given by fl9]). 
Only for P = (r_ = 0) we have Th = £-■ 



2.2 Hawking radiation of extremal black holes 

The extremal case can formally be considered as a limit of the non-extremal one, namely 
in the limit r + — > r_, i.e. Zq —¥ oo. In this limit we have 



cosh(7r(o;-A^)) cosh(7r((22;o-l)u;-A w /)) cosh(7r(w-A^)) e 



-2tt\ u 



z ->+oo cosh(7r(cu + A a; /)) cosh (7r( (2^0— l)w+Aa,/)) 



C0Sh(7T (CJ+Xui)) 



(33) 



Since one could doubt the legality of this limit because of the fact that r_ = r + is a singu- 
larity in our initial substitution f|T5|) for z, and also for completeness of our investigation, 
we will consider this case separately. 
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In the extremal case equation (I14p becomes 



d 
dr 



{r-r + ) Rui(r 



u 2 r 2 



(r — r+)'' 



1(1 + 1) 



ILi(r) = 0. 



(34) 



The smaller number of singular points is a significant reason to require a separate investi- 
gation, so we cannot expect to obtain again hypergeometric equation. In the extremal case 
our new variable will be: 



q:- 



r+ 



r — r + 



d_ 
dr 



r + dq 



(35) 



So the equation for Q w i(q) : = R w i(r) takes the form: 

q'QM + [u 2 -l(l + l) + 2u 2 q + u 2 q 2 ] Q wl (q) = 0. 



Here the same substitution \ u [ := a/w 2 — (1 + 1/2) 2 is appropriate for recognizing the 
above equation as Whittaker equation for W u i(z) := Qui(q), where z := 2iuq: 



W ul (z) + 



1/4 -{iK,iY , -^ 



+ 



W,Jz) = 0. 



(36) 



The general solution can be represented in terms of the confluent hypergeometric functions 
(Kummer functions [10J): 



W ut {z) = e - z/2 z^ +iX - 1 (d M(a, b, z) + C 2 U(a, b, z)) , where 



a = l/2 + i(w+A w j) 

b= l+2i\ ul . 



(37) 



The asymptotic solution on the horizon (r — y r + + =>- q — y +oo) follows the asymptotic 
expansion of the Kummer functions [10] : 



M(a,b,z)\ 



z — >+oo 



r(6) 



p ^zi7va 7 —a 



r(6-o) 



e * z a-6 



(o) n (H-0-fo) n 



IV. 



(-z)- n + 0(\z\- H ) 



Via) 



~R-1 

E 

_n=0 
5*— 1 
^ (6-a)n(l-a)n ( _ z) -n + (\ Z \- 8 ) 



+ 



n=0 



/?! 



where +ma is for — 7r/2<arg,2<37r/2 and — iira is for — 37r/2<arg,2< — tt/2. We also 
have 



U(a,b,z)\\ 



R-l 



J2 ia)n(1 ^, a ~ &)n (-^)-" + o(i^r^) 



n=0 



??■! 



3^ 



3^ 



y <arg^< — 



In our case the zero order terms in the sums (R = 1, S = 1) are sufficient again. Taking 
into account that argz = 7r/2 and also the expressions ( 1371) . one can obtain the solution on 
the horizon: 
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Q(g ->• +oo) ~ Ci r(lH-2iA 



Ull) 



le 



-f(A w i+w)„— iuiq 



(2i<jjqY 



+ 



e iG)q (2iuq) i 



T(l/2 - ifi-Ka)) T(l/2 + iiu+XJ)). 

+ C 2 e- i(bq (2iCjq) 



+ 



Here the suitable substitution that transforms the general solution in terms of ID wave 
modes in R is: 



2qe q 



-x/r 



-r Q r + r-r A 

x =- h r In ■ 



T — Ta 



2r A 



(38) 



which gives 



Rutir -> r+) ~ Ante*** + A^e 



A ou t '■— C\ 



T{l/2-i(u-\ ul )) 



+ C 2 e™ /2 c^, A m :=d 



r(l+2i\ u} i)e- s ?6u i 



r(l/2 + i(oj+\ ul ))- 
(39) 
For the asymptotic solution at the radial infinity (r — > +00 =>- q — > 0+ =>- z — > 0) we 
will use the relation between M(a, b, z) and [/(a, b, z) [TO] : 



£/(a, 6, z) 



Ti 



sin7r6 



M(a,b,z) 1 _ b M{(l+a-b,2-b,z) 



T(l+a-b)T(b) 



T(a)T(2-b) 



After applying the upper relation for z = where M(a, b, 0) = 1 we reach the asymptotic 
solution at the spatial infinity. 



Q(q ->. 0) ~ (2*wg) 



1/2 



Ci + - 



C<,i7r 



sinh(2vrA^)r(l/2 + i(u-\ ul ))T(l+2i\ ul ) 

C 2 m 



(2iojq) 



'A,, 



■(2iujq) 



-i\. 



sinh(27rA ajZ )r(l/2 + z(w+A wJ ))r(l-2tA w j) 
The substitutions (138]) and fc w j := A w //r again give the asymptotic solution in the 



form of ID wave in IR: 



R ul (r -> 00) 



J^± (B out e ik -' x + B m e~ ik ^) 



where 



Bout 
-Dill. ~ 



-C 2 m(2iCj) 1 ' 2 (iCj)- iX - 1 



sinh(27rA^)r(l/2 + i(u+\ ul ))T(l-2i\ wl ) 

C 2 iii 



C, + - 



sinh(2vrA^)r(l/2 + i(w-X ul ))T(l+2iX ul ) 



(40) 



(2iu) 



7M/2(iQ)i*ui. 



We should note that due to the time- reversal symmetry of wave equation it is permissible 
to work with the time-reversed formulation of the radiation boundary condition. In the 
normal picture the radiation boundary condition means missing of the in-going mode at the 
spatial infinity Bi n := 0. The reflection coefficient in this case is R := |^4m| 2 /|Antt| 2 - 

In the time-reversed picture the radiation mode becomes an in-going mode so in this 
treatment B out := 0. Respectively on the horizon the falling mode becomes an in-going 
mode and the reflected mode becomes out-going mode. So in this case we should take 

R := l-^-outl /I^T-ml • 

It is possible to turn out that 



A; 



A 



out 



+ 



A 



out 



Bi n :—0 



A, 



but always 



J^ou 



J\j 



A 



out 



Bi„:—0 



A 



out 



Aj r , 



£>ou 



In the work [9] the time-reversed picture is chosen despite the fact that there is no 
difference in the complexity of further calculations between both approaches. It is the same 
for our non-extremal black holes. But one could see from the expressions (1391) and (jlOj) that 
the time-reversed formulation gives a shorter way to R. All the results for R were made by 
us using both approaches for checking the correctness of all previous calculations. Following 
the time-reversed picture B out = =>- Ci = we find 



R 



\A 



out] 



I A I 2 

\- r Hn\ 



ie 


\ 7TUJ 


r(l+2«A w ,)w- 


-iui 




r(i/2 


-i{u-Xui)) 





r(i/2 + i(w+\,i)) 



2 0J % 



|r(i/2 + »(a+A Mt ))| a 2 « Kil 

~ |r(l/2-^(cu-A^))| 2 
For high frequencies Cj ^> / + 1/2 ( =^> X^i 



r(l+2zA^)e- 



cosh(7r(A^-a))) _ 2 ^ 



(41) 



-2iru> 



R 



oo 



-Attuj 



cosh(Tr(X LU i+u)) 
we have 



cosh(27rc<;) 



N zi e 



-i-KUJ 



(42) 



Hence we find the Hawking temperature in the extremal case in semiclassical approxi- 
mation 



T, 



H 



47rr 



(43) 



Contra-intuitively the temperature of the extremal case is non-zero and, as in the non- 
extremal case, is independent of the intrinsic characteristics of the black hole and depends 
only on the background parameter tq. 

3 Discussion 

In the present paper we studied the Hawking radiation of asymptotically non-flat dyonic 
black holes in 4D Einstein-Maxwell-dilaton gravity in semi-classical approximation. It was 
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shown that the problem can be solved exactly and we computed exactly the semi-classical 
radiation spectrum of both non-extremal and extremal black holes. 

Our results show that the Hawking temperature, calculated in the semi-classical approx- 
imation, is not compatible with the first law ffTUj) . The reason for this discrepancy is that the 
spacetime is asymptotically non-flat. Other examples for discrepancy between the Hawking 
temperature and the surface gravity in asymptotically non-flat spacetimes can be found 
in [TT], [12] and references therein. In principle, the relation between the Hawking radiation 
and the first law in asymptotically non-flat spacetimes is controversial and depends on the 
particular case. In general, the first law in asymptotically non-flat spacetime is not directly 
connected to the temperature of the particle flux at infinity. For example, the black holes 
considered in the present work can not emit massive particles because the mass term in the 
wave equation ( TTTj) leads to the appearance of confining potential (growing unboundedly to 
infinity) which prevents the particles from escaping to infinity. In other words, the Hawking 
temperature for massive particles, measured by an asymptotic observer, is zero. 

The asymptotical non-flatness leads to ambiguous thermodynamical characteristics and 
thermodynamics as a whole. For example, the surface gravity is given by the formula 

ev M r = «r (44) 

on the horizon. The above definition, however, gives the surface gravity up to a constant 
because there is a freedom to rescale £ by a constant. In the asymptotically flat case this 
rescaling freedom can be fixed by choosing a unit norm for the Killing field at infinity. In 
linear dilaton spacetimes there is no natural way to fix the rescaling freedom. The choice 
£ = d/dt made in the papers devoted to the linear dilaton black holes is thus ad hoc. The 
ambiguity in choosing the time vector field is present also in the Hamiltonian formalism for 
linear dilaton spacetimes [5] where the time vector field is again £ = d/dt. 

One possible way to overcome the problem of the rescaling freedom in our case is to 
replace the Killing field £ with the Kodama vector field K = v/y^ which has a unit norm 
at infinity [13]. This investigation is in progress and the results will be presented elsewhere. 
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